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We report first-principles calculations of the frequency dependent linear and second-order optical 
properties of the two probable extended-solid phases of CO2-V, i.e. 742d and P2i2i2i. Compared 
to the parent Cmca phase the linear optical susceptibility of both phases is much smaller. We 
find that I42d and P2i2i2i differ substantially in their linear optical response in the higher energy 
regime. The nonlinear optical responses of the two possible crystal structures differ by roughly a 
factor of five. Since the differences in the nonlinear optical spectra are pronounced in the low energy 
regime, i.e. below the band gap of diamond, measurements with the sample inside the diamond 
anvil cell are feasible. We therefore suggest optical experiments in comparison with our calculated 
data as a tool for the unambiguous identification of the high pressure phase of CO2. 

PACS numbers: 61.50A, 42.65K, 61.50K, 78.66H 



A few years after the synthesis of the extended-solid 
phase of CO2 at 40 GPa the crystal structure of this 
phase is still not known. This phase was found to be 
similar to one of the noncentrosymmetric polymorphs of 
quartz 0, Q and was designated as CO2-V. Since solid 
CO2 has very high tensile strength (comparable to that 
of boron nitride y|) this discovery could be very useful 
in device applications. There is also considerable interest 
in the role of CO2 in planetary geology p|. 

Experimental work done to identify the exact crystal 
structure of CO2 at high pressure suggests a low sym- 
metry tridymite phase (P2i2i2i) [5J. The presence of a 
nonlinear optical signal also indicated the absence of in- 
version symmetry in the crystal structure . Most of the 
theoretical work has been concentrated around the lat- 
tice relaxations and total energy calculations in search 
of the minimum energy structure [(J, 0, @ ■ Theory pre- 
dicts the /3-cristobalite phase (742c?) to be lower in energy 
than P2i2i2i by around 0.19 eV per formula unit Q. 
However, none of the Si02-like polymorphs of COoKave 
a satisfactory fit to the experimental XRD data ja, @ • 
Other possible quartz polymorphs have been eliminated 
by their much higher total energies 0. Since all experi- 
ments were performed at elevated temperatures (1800K) 
and the theoretical energy difference between the two 
phases is small, the energy minimization is not sufficient 
to solve the discrepancy between theory and experiments 
and to predict the actual high pressure phase. In order 
to understand the nature of this phase transition there 
is the need for analyzing other physical properties and 
finding a satisfactory agreement between theory and ex- 
periment. Since the two probable phases have very differ- 
ent electronic structure it is possible to distinguish them 
by their optical response. So far, neither a theoretical 
or experimental study of the linear optical response have 
been performed. Although second harmonic generation 
(SHG) by CO2 has been seen experimentally, there exist 
no experimental data or calculations of the frequency de- 



pendent SHG susceptibility. The aim of the present work 
is to show that the linear and nonlinear optical spectra 
indeed allow an unambiguous identification of the CO2- 
V crystal structure. To this extent we have calculated the 
linear and second-order optical response of solid carbon 
dioxide, CO2-V, in the two possible phases. The linear 
optical spectra of these probable phases are compared to 
each other and also to that of the Cmca phase, which 
under high pressure and temperature conditions leads 
to CO2-V. Since quartz is often used as a standard in 
second-harmonic (SH) experiments, a detailed compari- 
son is made between the second-order optical response of 
the two CO2-V phases and that of SiC>2 (IA2d phase). 
We propose that on basis of our results, future exper- 
imental work will help to resolve the present disagree- 
ment between theory and experiment and hence clarify 
the high pressure phase of CO2. 

Total energy calculations are performed using the 
state-of-the-art full-potential linearized augmented plane 
wave FPLAPW method implemented in the WIEN2k 
code 0. For the exchange-correlation potential we use 
the generalized gradient approximation (GGA) derived 
by Perdew and Wang [lOj . The detailed formalism for the 
determination of the linear susceptibility tensor \ 
within the FPLAPW formalism has been presented be- 
fore [HI. The susceptibility for the SHG x (2) ( 2 ^ w , w ) 
has been calculated using an extension to this program 
(l2^ . The SHG susceptibility obtained satisfies all the 
theoretical sum rules presented by Scandolo and Bassani 
[l3| . All the calculations are converged in terms of basis 
functions as well as in the size of the fc-point mesh rep- 
resenting the Brillouin zone. The optical properties are 
calculated on a mesh of 1500 fc-points in the irreducible 
wedge of the Brillouin zone. 

The phases of CO2 under investigation in the present 
work are the Cmca phase and the two probable extended- 
solid phases for CO2-V, namely I42d and P2i2i2i. For 
the P2i2i2i phase the starting point is the experimen- 
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tal lattice parameters |5( with the atomic positions of 
the C222i phase (which is the supergroup of the space 
group P2 1 2 1 2 1 ). In a first step, the lattice parameters 
are kept constant at the experimental values and the in- 
ternal parameters are fully relaxed. This structure is 
referred to as P2i2i2i~expt in the present work. Next, 
the lattice and the internal parameters are simultane- 
ously relaxed (keeping a constant volume Q) to attain 
the minimum energy configuration. This structure is re- 
ferred to as P2\2\l\-rel. The energy lowering on going 
from P2\2\2\-ex'pt to the P2\2\l\-rel structure is 0.81 
eV per formula unit, which is around 0.16 eV more than 
the previous theoretical value Q. The larger energy in 
the present calculations compared to the previous theo- 
retical value could be due to the difference in the method 
of calculation, since in theprevious work the pseudopo- 
tential approach was used y|. In case of the IA2d phase, 
due to the lack of experimental data, the starting lat- 
tice constants and internal parameters are taken to be 
those of SiC>2 in the same phase. Then the atoms as 
well as the cell volume are fully relaxed to attain the 
minimum energy configuration. The lattice parameters 
for all the structures are presented in Table I. In accor- 
dance with the previous theoretical work 0, Q we find 
/?-cristobalite (IA2d) to be lower in energy compared to 
the P2i2i2i~rel phase. The difference in energy between 
the P2i2i2i~rel and the theoretically determined I42d 
phase at 40 GPa is 0.38 eV/C0 2 . 
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FIG. 1: The magnitude of the linear susceptibility compo- 
nents for solid CO2 in (a) the Cmca phase, (b) the I42d 
phase at 40 GPa, (c) the P2i2i2i-expt structure and (d) the 
P2i2i2i-rei structure. 

The results for the linear optical spectra of CO2 in 



TABLE I: Lattice parameters for various phases. The cell 
parameters are in A. The atomic positions are in lattice co- 
ordinates. 



Phase 


a (A) b (A) c (A) 


atom position (a,b,c) 


J42d 


4.162 4.162 5.886 


C (0,0,0) 

O (0.880,0.250,0.125) 


P2 1 2 1 2 1 -expt 


6.214 4.351 6.064 


C(l) (0.434,0.620,0.690) 
C(2) (0.920,0.940,0.710) 
0(1) (0.890,0.599,0.200) 
0(2) (0.550,0.200,0.020) 
0(3) (0.228,0.577,0.214) 
0(4) (0.469,0.257,0.351) 


P2 1 2 1 2 1 -rel 


6.767 3.773 6.420 


C(l) (0.900,0.920,0.700) 
C(2) (0.420,0.494,0.690) 
0(1) (0.589,0.190,0.990) 
0(2) (0.247,0.703,0.211) 
0(3) (0.430,0.220,0.293) 
0(4) (0.934,0.757,0.161) 



various phases are presented in Fig. ^ None of the opti- 
cal spectra are scissors corrected. Both the P2\2\2\ and 
the l\2d phase are taken at 40 GPa. The Cmca phase 
shows a strong anisotropy with \Xyy ' ( w )| much larger in 
amplitude than |Xxx(w)|, and |xii ( w )| lying in between 
the two. The reason for this anisotropy is the orientation 
of the molecules within the crystal. The neighbouring 
molecules are parallel to each other in the ab plane and 
form a herringbone pattern in the ac and be planes [Fig. 
2(b) in Ref. [l5j. This causes the yy and zz polarized 
transition matrix elements (TMEs) to be large. On the 
other hand the xx polarized TMEs are close to zero, lead- 
ing to small values for |xix (k>)|- \xJz ■ (w)| is smaller than 
\Xm (k>)l because the lattice parameter b is smaller than 
c. On going from the low pressure Cmca phase to any 
other phase of CO2-V, \\y~y (w) and the anisotropy be- 
tween (w)| and Ixiz ( w )l decrease dramatically. The 
latter fact can again be explained by the orientation of 
the molecules inside the unit cell [Fig. 3E and 3F in Ref. 
8J. The neighbouring CO2 molecules are no more parallel 
or perpendicular in any direction leading to a weakening 
of the very strong transitions and an enhancement of the 
TMEs which were forbidden in the Cmca phase. For 
all the probable phases of CO2-V the spectral weights 
are shifted to higher energies. The optical band gap of 
the P2i2i2i~expt structure is smaller than that of the 
Cmca phase, while the other two structures show almost 
equal optical band gaps, which are higher than that of 
the Cmca phase. Compared to the Cmca phase the max- 
imum peak height is around 3 times smaller in the /42<i 
phase, while it is about 4 times smaller in the P2i2i2i 
structures. Beyond 11 eV, the response of all the proba- 
ble phases of CO2-V is on average similar in magnitude, 
but differs in the amount of anisotropy. 



For I42d there are two independent tensor compo- 
nents. Compared to the zz component, the xx compo- 
nent is approximately upto 2 times as large in the energy 
range below 10 eV. Upto 10.5 eV and between 11.5 and 
13 eV the effect is opposite, while the two are nearly 
the same for higher energies. For the P2\2\2\ struc- 
tures there are three incquivalcnt components and the 
anisotropy between the xx and zz components is very 
small compared to the I42d phase. We also find the lat- 
tice relaxation effects which have already been described 
in detail by Dong et al. 0. They have some conse- 
quences for the optical spectra. On going from the Cmca 
phase to P2 1 2 1 2 1 -expt, the optical band gap decreases, 
while it increases in case of P2 1 2 1 2 1 -reL Another notice- 
able feature is that the anisotropy is much larger in the 
P2 1 2i2i~expt structure. Since CO2-V is not quenchable 
into ambient pressure at room temperature the optical 
measurements need to be performed with the sample in- 
side the diamond anvil cell (DAC). The band gap of dia- 
mond is 5.5eV. Therefore the energy regime below 5.5eV 
is most interesting, since the optical response in this en- 
ergy regime is purely from the sample inside the DAC. 
The difference between the linear optical spectra of vari- 
ous probable phases of CO2 below this energy range is not 
as pronounced as in the higher energy regime. The results 
though are of theoretical interest and help in identifica- 
tion of various features in the nonlinear optical spectra 
(which is rather difficult to do from the band structure 
plots), but may not be very good in phase distinction 
under these experimental conditions. 
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FIG. 2: The magnitude of the second-harmonic generation 
susceptibility components xyz and zyx for (a) Si02 in the 
J42d phase, (b) CO2 in the 742d phase at 40 GPa, (c) in the 
P2\2\2\-expt structure and (d) in the P2i2i2i-reZ structure. 



The nonlinear optical properties are much more sen- 
sitive to small changes in the band structure than the 
linear optical spectra. The calculated magnitude of the 
SHG susceptibility by CO2-V in various probable phases 
are presented in Fig. |2 Due to the presence of inversion 
symmetry in the Cmca phase the SHG susceptibility is 
zero. The second-order response of Si02 (I42d), which is 
often used as standard in SHG experiments, is also pre- 
sented in the figure. In the energy range between 3 and 
5.5 eV the magnitude of the SHG susceptibility by CO2 in 
the P2\2\2\-expt phase is 1.5-2 times smaller than that 
of quartz. On the other hand the IA2d phase of CO2 
shows a SHG susceptibility 3-4 times larger than that of 
Si02- These pronounced differences make the phases eas- 
ily identifiable. The energy range of 3-5.5 eV lies below 
the band gap of diamond and corresponds to laser light 
of wavelengths between 413.3 - 247.9 nm. Previous SH 
experiments were performed using laser light of 1054 nm 
wavelength 0, which correspond to 1.176 eV. At this 

energy our calculated value of \x\ (2w,w,w)| for CO2 
in the I42d phase is 2.7 times larger than that of Si02- 
While, the SHG susceptibility by C0 2 in the P2 1 2 x 2 1 - 
expt phase is 2.6 times smaller than that of Si02- The 
previous experimental work [j| reports the SH conversion 
efficiency of the CO2-V phase, which is not just propor- 
tional to the square of the SHG susceptibility but also 
depends on the input intensity and sample parameters. 
So no direct comparison of the existing experimental data 
can be made with the present results. Since the differ- 
ence in the second-order response by the probable phases 
of CO2-V are more pronounced in the energy range be- 
tween 3 and 5.5 eV than in between 1-2 eV, experiments 
at higher energies would be favourable. Experimental 
determination of the SHG susceptibility by CO2-V and 
Si02, over some frequency range, and a comparison with 
our work could lead to conclusive results. 

The two structures of the P2\2\2\ phase are also differ- 
entiable by their SHG response. Due to the large optical 
band gap of P2\2\2\- J rel the SHG susceptibility is ap- 
proximately 3 times smaller than that of P2\2\2\-expt. 
As shown in Ref. 7 the XRD patterns for neither of 
these structures matches the experiments. The theoret- 
ically calculated XRD pattern of P2\2\2\-rel shows a 
better agreement with the experimental XRD pattern at 
higher values of 20 (9-13 degrees), while the XRD pattern 
for P2\2\2\-expt is in better agreement at low values of 
29 (6-8 degrees). It is possible that under the experi- 
mental conditions of high temperatures the material is 
a polycrystalline mixture of these two structures or has 
a configuration in between the two, in which case the 
difference between the optical spectra of the two phases 
would represent a measure of the uncertainty in the op- 
tical response of this phase. Despite this fact, we want 
to point out that in general the differences in the op- 
tical response between the two phases are much more 
pronounced than "the lattice relaxation effects" for the 
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P2i2i2i phase making them distinguishable by optical 
probes. 

The optical coefficients in the static limit are directly 
measurable quantities. The calculated inequivalent com- 
ponents of ei(0) and Ix*- 2 - 1 (0) | are presented in the Table 
II. 

TABLE II: The static dielectric constants ef (0), ef (0) and 
ef z (0), and the SHG susceptibility components in the static 

limit |xii/^(0)| and |xS/z(0)l m xiO -8 esu for CO2 in various 
phases. 



Phase 


ef (0) ef (0) ef (0) 


lxil(o)| |xSUo)| 


Cmca 


1.87 3.23 2.84 


0.00 0.00 


U2d 


2.82 2.82 2.66 


0.74 0.69 


P2 1 2 1 2i-expt 


2.69 3.57 2.67 


0.10 0.25 


P2 1 2 1 2 1 -rel 


2.69 2.86 2.68 


0.02 0.04 



In the Cmca phase, the anisotropy between ef (0) and 
ef (0) is_54%, and between ef (0) and ef (0) it is 41%. 
In the 742d phase, the anisotropy between the xx and zz 
components is 6% only. For the tridymite phase at exper- 
imental lattice parameters (P2 1 2 1 2 1 -ea;pt) the maximum 
anisotropy found between ef (0) and ef (0) (28 %), while 
between ef (0) and ef (0) it is only 1%. In the case of 
P2i2i2i-re/ the anisotropy between the xx and yy com- 
ponents is 7%, and between the xx and zz components 
is close to zero. 

Amongst the phases under investigation in the present 

work the anisotropy in the xyz and zyx components of 

|X^(0)l is maximum for P2\2i2i-expt with 86 % and 

minimum for 742<i with 7%. P2\2\2\-rel lies in between 

the two with an anisotropy of 67 %. For SiC>2 in the I42d 

(2) (2) 
phase, our calculated values of |Xxj/z(0)l an d IXzjMO)! 

are 0.28 and 0.20 xlO -8 esu, respectively. It may be 

noted that the anisotropy as well as the magnitude of 

(2) 

Xi (0) for CO2 in the P2 1 2 1 2i~expt phase are close to 
that of Si0 2 . The fact that |x (2) (0)| of C0 2 is around 3 
times larger in the 742c? phase, and 5-10 times smaller 
in the P2\2i2i-rel phase compared to Si02 makes the 
distinction between the two phases clear and hence the 
identification of the high pressure phase feasible. 

In summary we find that the optical properties are ex- 
cellent tools for an unambiguous determination of the 
high pressure phase of CO2. Particularly, the second- 
order response of CO2-V compared to that of SiC>2 can be 
used for phase determination even with the sample inside 



the DAC. Experiments of this kind and a detailed com- 
parison with the present theoretical results would shed 
light on the existing discrepancy between theory and ex- 
periments regarding the space group of the CO2-V phase. 
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We report first-principles calculations of the frequency dependent linear and second-order optical 
properties of the two probable extended-solid phases of CO2-V, i.e. IA2d and P2i2i2i. Compared 
to the parent Cmca phase the linear optical susceptibility of both phases is much smaller. We 
find that IA2d and P2i2i2i differ substantially in their linear optical response in the higher energy 
regime. The nonlinear optical responses of the two possible crystal structures differ by roughly a 
factor of five. Since the differences in the nonlinear optical spectra are pronounced in the low energy 
regime, i.e. below the band gap of diamond, measurements with the sample inside the diamond 
anvil cell are feasible. We therefore suggest optical experiments in comparison with our calculated 
data as a tool for the unambiguous identification of the high pressure phase of CO2 . 

PACS numbers: 61.50A, 42.65K, 61.50K, 78.66H 



A few years after the synthesis of the extended-solid 
phase of CO2 at 40 GPa the crystal structure of this 
phase is still not known. This phase was found to be 
similar to one of the noncentrosymmetric polymorphs of 
quartz [1, 2] and was designated as CO2-V. Since solid 
CO2 has very high tensile strength (comparable to that 
of boron nitride [3]) this discovery could be very useful 
in device applications. There is also considerable interest 
in the role of CO2 in planetary geology [4]. 

Experimental work done to identify the exact crystal 
structure of CO2 at high pressure suggests a low sym- 
metry tridymite phase (P2i2i2i) [5]. The presence of a 
nonlinear optical signal also indicated the absence of in- 
version symmetry in the crystal structure [1]. Most of the 
theoretical work has been concentrated around the lat- 
tice relaxations and total energy calculations in search of 
the minimum energy structure [6-8]. Theory predicts 
the /3-cristobalite phase (I42d) to be lower in energy 
than P2\l\l\ by around 0.19 eV per formula unit [7]. 
However, none of the SiCVlike polymorphs of CO2 gave 
a satisfactory fit to the experimental XRD data [5, 8]. 
Other possible quartz polymorphs have been eliminated 
by their much higher total energies [6]. Since all experi- 
ments were performed at elevated temperatures (1800K) 
and the theoretical energy difference between the two 
phases is small, the energy minimization is not sufficient 
to solve the discrepancy between theory and experiments 
and to predict the actual high pressure phase. In order 
to understand the nature of this phase transition there 
is the need for analyzing other physical properties and 
finding a satisfactory agreement between theory and ex- 
periment. Since the two probable phases have very differ- 
ent electronic structure it is possible to distinguish them 
by their optical response. So far, neither a theoretical 
or experimental study of the linear optical response have 
been performed. Although second harmonic generation 
(SHG) by CO2 has been seen experimentally, there exist 
no experimental data or calculations of the frequency de- 



pendent SHG susceptibility. The aim of the present work 
is to show that the linear and nonlinear optical spectra 
indeed allow an unambiguous identification of the CO2- 
V crystal structure. To this extent we have calculated the 
linear and second-order optical response of solid carbon 
dioxide, CO2-V, in the two possible phases. The linear 
optical spectra of these probable phases are compared to 
each other and also to that of the Cmca phase, which 
under high pressure and temperature conditions leads 
to CO2-V. Since quartz is often used as a standard in 
second-harmonic (SH) experiments, a detailed compari- 
son is made between the second-order optical response of 
the two CO2-V phases and that of SiC>2 (I42d phase). 
We propose that on basis of our results, future exper- 
imental work will help to resolve the present disagree- 
ment between theory and experiment and hence clarify 
the high pressure phase of CO2. 

Total energy calculations are performed using the 
state-of-the-art full-potential linearized augmented plane 
wave FPLAPW method implemented in the WIEN2k 
code [9]. For the exchange-correlation potential we use 
the generalized gradient approximation (GGA) derived 
by Perdew and Wang [10]. The detailed formalism for the 
determination of the linear susceptibility tensor %W (w) 
within the FPLAPW formalism has been presented be- 
fore [11]. The susceptibility for the SHG x (2) ( 2w > w > w ) 
has been calculated using an extension to this program 
[12]. The SHG susceptibility obtained satisfies all the 
theoretical sum rules presented by Scandolo and Bassani 
[13]. All the calculations are converged in terms of basis 
functions as well as in the size of the fc-point mesh rep- 
resenting the Brillouin zone. The optical properties are 
calculated on a mesh of 1500 fc-points in the irreducible 
wedge of the Brillouin zone. 

The phases of CO2 under investigation in the present 
work are the Cmca phase and the two probable extended- 
solid phases for C0 2 -V, namely I42d and P2i2i2i. For 
the _P2i2i2i phase the starting point is the experimen- 



tal lattice parameters [5] with the atomic positions of 
the C222i phase (which is the supergroup of the space 
group P2i2i2i). In a first step, the lattice parameters 
are kept constant at the experimental values and the in- 
ternal parameters are fully relaxed. This structure is 
referred to as P2i2i2i-expt in the present work. Next, 
the lattice and the internal parameters are simultane- 
ously relaxed (keeping a constant volume [14]) to attain 
the minimum energy configuration. This structure is re- 
ferred to as P2i2i2i-reL The energy lowering on going 
from P2 1 2 1 2 1 -expt to the P2 1 2 1 2 1 -rel structure is 0.81 
eV per formula unit, which is around 0.16 eV more than 
the previous theoretical value [7]. The larger energy in 
the present calculations compared to the previous theo- 
retical value could be due to the difference in the method 
of calculation, since in the previous work the pseudopo- 
tential approach was used [7]. In case of the I42d phase, 
due to the lack of experimental data, the starting lat- 
tice constants and internal parameters are taken to be 
those of SiC>2 in the same phase. Then the atoms as 
well as the cell volume are fully relaxed to attain the 
minimum energy configuration. The lattice parameters 
for all the structures are presented in Table I. In accor- 
dance with the previous theoretical work [7, 8] we find 
/?-cristobalite (I42d) to be lower in energy compared to 
the P2\2\2\-rel phase. The difference in energy between 
the P2\2\2\-rel and the theoretically determined 142d 
phase at 40 GPa is 0.38 eV/C0 2 . 
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FIG. 1: The magnitude of the linear susceptibility compo- 
nents for solid CO2 in (a) the Cmca phase, (b) the I42d 
phase at 40 GPa, (c) the P2i2i2\-expt structure and (d) the 
P2i2i2i-re/ structure. 

The results for the linear optical spectra of CO2 in 



TABLE I: Lattice parameters for various phases. The cell 
parameters are in A. The atomic positions are in lattice co- 
ordinates. 



Phase 


a (A) 6 (A) c (A) 


atom position (a,b,c) 


J42d 


4.162 4.162 5.886 


C (0,0,0) 

O (0.880,0.250,0.125) 


P2 1 2 1 2 1 -expt 


6.214 4.351 6.064 


C(l) (0.434,0.620,0.690) 
C(2) (0.920,0.940,0.710) 
0(1) (0.890,0.599,0.200) 
0(2) (0.550,0.200,0.020) 
0(3) (0.228,0.577,0.214) 
0(4) (0.469,0.257,0.351) 


P2 1 2 1 2 1 -rel 


6.767 3.773 6.420 


C(l) (0.900,0.920,0.700) 
C(2) (0.420,0.494,0.690) 
0(1) (0.589,0.190,0.990) 
0(2) (0.247,0.703,0.211) 
0(3) (0.430,0.220,0.293) 
0(4) (0.934,0.757,0.161) 



various phases are presented in Fig. 1. None of the opti- 
cal spectra are scissors corrected. Both the P2\2\2\ and 
the 142d phase are taken at 40 GPa. The Cmca phase 
shows a strong anisotropy with Ixs/V ( w )l much larger in 
amplitude than and |xi^( w )l lying m between 

the two. The reason for this anisotropy is the orientation 
of the molecules within the crystal. The neighbouring 
molecules are parallel to each other in the ab plane and 
form a herringbone pattern in the ac and be planes [Fig. 
2(b) in Ref. 15]. This causes the yy and zz polarized 
transition matrix elements (TMEs) to be large. On the 
other hand the xx polarized TMEs are close to zero, lead- 
ing to small values for |xia;(w)|. |Xzz ( w )l l& smaller than 
\Xyy (pS)\ because the lattice parameter b is smaller than 
c. On going from the low pressure Cmca phase to any 
other phase of CO2-V, \xyy (w)| and the anisotropy be- 
tween \x^yy (w)| and |xia;(w)| decrease dramatically. The 
latter fact can again be explained by the orientation of 
the molecules inside the unit cell [Fig. 3E and 3F in Ref. 
8]. The neighbouring CO2 molecules are no more parallel 
or perpendicular in any direction leading to a weakening 
of the very strong transitions and an enhancement of the 
TMEs which were forbidden in the Cmca phase. For 
all the probable phases of CO2-V the spectral weights 
are shifted to higher energies. The optical band gap of 
the P2i2i2i-expt structure is smaller than that of the 
Cmca phase, while the other two structures show almost 
equal optical band gaps, which are higher than that of 
the Cmca phase. Compared to the Cmca phase the max- 
imum peak height is around 3 times smaller in the I42d 
phase, while it is about 4 times smaller in the P2\2]2\ 
structures. Beyond 11 eV, the response of all the proba- 
ble phases of CO2-V is on average similar in magnitude, 
but differs in the amount of anisotropy. 



For I42d there are two independent tensor compo- 
nents. Compared to the zz component, the xx compo- 
nent is approximately upto 2 times as large in the energy 
range below 10 eV. Upto 10.5 eV and between 11.5 and 
13 eV the effect is opposite, while the two are nearly 
the same for higher energies. For the P2i2i2i struc- 
tures there are three inequivalent components and the 
anisotropy between the xx and zz components is very 
small compared to the I42d phase. We also find the lat- 
tice relaxation effects which have already been described 
in detail by Dong et al. [7]. They have some conse- 
quences for the optical spectra. On going from the Cmca 
phase to P2i2i2i-expt, the optical band gap decreases, 
while it increases in case of P2\2\2\-rel. Another notice- 
able feature is that the anisotropy is much larger in the 
P2\2\2\-exTpt structure. Since CO2-V is not quenchable 
into ambient pressure at room temperature the optical 
measurements need to be performed with the sample in- 
side the diamond anvil cell (DAC). The band gap of dia- 
mond is 5.5eV. Therefore the energy regime below 5.5eV 
is most interesting, since the optical response in this en- 
ergy regime is purely from the sample inside the DAC. 
The difference between the linear optical spectra of vari- 
ous probable phases of CO2 below this energy range is not 
as pronounced as in the higher energy regime. The results 
though are of theoretical interest and help in identifica- 
tion of various features in the nonlinear optical spectra 
(which is rather difficult to do from the band structure 
plots), but may not be very good in phase distinction 
under these experimental conditions. 
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FIG. 2: The magnitude of the second-harmonic generation 
susceptibility components xyz and zyx for (a) SiCh in the 
I42d phase, (b) C0 2 in the J42d phase at 40 GPa, (c) in the 
P2\2i2i-expt structure and (d) in the P2i2i2i-re/ structure. 



The nonlinear optical properties are much more sen- 
sitive to small changes in the band structure than the 
linear optical spectra. The calculated magnitude of the 
SHG susceptibility by CO2-V in various probable phases 
are presented in Fig. 2. Due to the presence of inversion 
symmetry in the Cmca phase the SHG susceptibility is 
zero. The second-order response of SiC"2 (I42d), which is 
often used as standard in SHG experiments, is also pre- 
sented in the figure. In the energy range between 3 and 
5.5 eV the magnitude of the SHG susceptibility by CO2 in 
the P2\2i2i-expt phase is 1.5-2 times smaller than that 
of quartz. On the other hand the IA2d phase of CO2 
shows a SHG susceptibility 3-4 times larger than that of 
Si02- These pronounced differences make the phases eas- 
ily identifiable. The energy range of 3-5.5 eV lies below 
the band gap of diamond and corresponds to laser light 
of wavelengths between 413.3 - 247.9 nm. Previous SH 
experiments were performed using laser light of 1054 nm 
wavelength [1], which correspond to 1.176 eV. At this 
energy our calculated value of |xi 2 '(2 )\ for C0 2 

in the I42d phase is 2.7 times larger than that of Si02- 
While, the SHG susceptibility by C0 2 in the P2i2i2i- 
expt phase is 2.6 times smaller than that of Si02- The 
previous experimental work [5] reports the SH conversion 
efficiency of the CO2-V phase, which is not just propor- 
tional to the square of the SHG susceptibility but also 
depends on the input intensity and sample parameters. 
So no direct comparison of the existing experimental data 
can be made with the present results. Since the differ- 
ence in the second-order response by the probable phases 
of CO2-V are more pronounced in the energy range be- 
tween 3 and 5.5 eV than in between 1-2 eV, experiments 
at higher energies would be favourable. Experimental 
determination of the SHG susceptibility by CO2-V and 
Si02, over some frequency range, and a comparison with 
our work could lead to conclusive results. 

The two structures of the P2\2\2\ phase are also differ- 
entiable by their SHG response. Due to the large optical 
band gap of P2\2\2\-rel the SHG susceptibility is ap- 
proximately 3 times smaller than that of P2\2\2\-expt. 
As shown in Ref. 7 the XRD patterns for neither of 
these structures matches the experiments. The theoret- 
ically calculated XRD pattern of P2\2\2\-rel shows a 
better agreement with the experimental XRD pattern at 
higher values of 29 (9-13 degrees), while the XRD pattern 
for P2\2\2\-expt is in better agreement at low values of 
29 (6-8 degrees). It is possible that under the experi- 
mental conditions of high temperatures the material is 
a polycrystalline mixture of these two structures or has 
a configuration in between the two, in which case the 
difference between the optical spectra of the two phases 
would represent a measure of the uncertainty in the op- 
tical response of this phase. Despite this fact, we want 
to point out that in general the differences in the op- 
tical response between the two phases are much more 
pronounced than "the lattice relaxation effects" for the 



P2\l\l\ phase making them distinguishable by optical 
probes. 

The optical coefficients in the static limit are directly 
measurable quantities. The calculated inequivalent com- 
ponents of ei(0) and |x' 2 '(0)l are presented in the Table 
II. 

TABLE II: The static dielectric constants ef x (0), e\ v [fi) and 
e| z (0), and the SHG susceptibility components in the static 
limit |xz»a;(0)| and |xa:»z(0)| in xlO esu for CO2 in various 
phases. 



Phase 


ef"(0) ef(0) eV(0) 


1x^(0)1 lx^(o)| 


Cmca 


1.87 3.23 2.84 


0.00 0.00 


IA2d 


2.82 2.82 2.66 


0.74 0.69 


P2 1 2 1 2 1 -expt 


2.69 3.57 2.67 


0.10 0.25 


P2 1 2 1 2 1 -rel 


2.69 2.86 2.68 


0.02 0.04 



In the Cmca phase, the anisotropy between ef x (0) and 
ef (0) is_54%, and between ef x (0) and ef (0) it is 41%. 
In the I42d phase, the anisotropy between the xx and zz 
components is 6% only. For the tridymite phase at exper- 
imental lattice parameters (P2\2\2\-expt) the maximum 
anisotropy found between ef v (0) and ef z (0) (28 %), while 
between ef x (0) and ef z (0) it is only 1%. In the case of 
P2\2\2\-rel the anisotropy between the xx and yy com- 
ponents is 7%, and between the xx and zz components 
is close to zero. 

Amongst the phases under investigation in the present 
work the anisotropy in the xyz and zyx components of 
|x' 2 '(0)l m maximum for P2i2i2i-expt with 86 % and 
minimum for I42d with 7%. P2\2\2\-rel lies in between 

the two with an anisotropy of 67 %. For SiC-2 in the I42d 

(2) (2) 
phase, our calculated values of |xxj/z(0)| and |xzyai(0)| 

are 0.28 and 0.20 xl0~ 8 esu, respectively. It may be 
noted that the anisotropy as well as the magnitude of 
x\ (0) for CO2 in the P2i2i2i-expt phase are close to 
that of Si0 2 . The fact that |x (2) (0)| of C0 2 is around 3 
times larger in the I42d phase, and 5-10 times smaller 
in the P2i2i2i-re/ phase compared to SiC-2 makes the 
distinction between the two phases clear and hence the 
identification of the high pressure phase feasible. 

In summary we find that the optical properties are ex- 
cellent tools for an unambiguous determination of the 
high pressure phase of CO2. Particularly, the second- 
order response of CO2-V compared to that of SiC-2 can be 
used for phase determination even with the sample inside 



the DAC. Experiments of this kind and a detailed com- 
parison with the present theoretical results would shed 
light on the existing discrepancy between theory and ex- 
periments regarding the space group of the CO2-V phase. 
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